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We study potentially observable consequences of spatiotemporal discreteness for the motion of 
massive and massless particles. First we describe some simple models for the motion of a massive 
point particle in a fixed causal set background. If the causal set is faithfully embeddable in Minkoswki 
spacetime, the models give rise to particle motion in the continuum spacetime. At large scales, the 
microscopic swerves induced by the underlying atomicity manifest themselves as a Lorentz invariant 
diffusion in energy-momentum governed by a single phenomenological parameter, and we derive in 
full the corresponding diffusion equation. Inspired by the simplicity of the result, we then derive the 
most general Lorentz invariant diffusion equation for a massless particle, which turns out to contain 
two phenomenological parameters describing, respectively, diffusion and drift in the particle's energy. 
The particles do not leave the light cone however: their worldlines continue to be null geodesies. 
Finally, we deduce bounds on the drift and diffusion constants for photons from the blackbody 
nature of the spectrum of the cosmic microwave background radiation. 

PACS numbers: 04.60.Bc 



I. INTRODUCTION 

The search for a theory of quantum gravity is not, as 
yet, motivated by experimental results. We currently 
have no unambiguously relevant quantum gravitational 
phenomena to guide us in developing candidate theo- 
ries, though it has long been suggested that a nonzero 
value of the cosmological constant of order 10~^^° could 
have a quantum gravitational origin [l|, [3, H, 01 • Outside 
of cosmology, black hole thermodynamics is often men- 
tioned as one example of a realm where concepts of gen- 
eral relativity and quantum mechanics must both come 
into play - but experimental black hole physics is out of 
our reach for now and even if analogue models of black 
holes in condensed matter systems could be tested this 
would only probe the semiclassical regime and not full 
quantum gravity. The existing approaches to quantum 
gravity have therefore been developed with the hope that 
the confrontation with experiment can be postponed. At 
the present time, however, the growing number of dif- 
ferent approaches means that the importance of testing 
ideas against observation, if at all possible, is greater than 
ever. 

Experimental verification of quantitative and unex- 
pected predictions is of the utmost importance in the 
development of a successful new theory. An example per- 
tinent to the current paper is if we were to find observa- 
tional evidence that spacetime is fundamentally discrete, 
then that would have a major impact on the direction of 
quantum gravity research. What form might such evi- 
dence take; what could be the Brownian motion of our 
age? To answer that question requires the development 
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of phenomenology that draws on essential aspects of a 
discrete theory of quantum gravity which turns out to be 
achievable in the causal set approach. 

Causal set theory is a discrete, Lorentz invariant ap- 
proach to quantum gravity [1, [a, ■ For reviews and 
further references see, for example, K is a work 

in progress: a quantum causal set dynamics still eludes 
us. Without a quantum dynamics it seems at first sight 
premature to develop causal set phenomenology but the 
kinematics of causal set theory is so concrete that we are 
able to make some progress in this direction. 

A causal set is a locally finite partial order and is the 
kinematical basis for the theory. One could state the 
central hypothesis as that spacetime is a causal set, or, if 
one wanted to hedge one's bets whilst the foundations of 
quantum theory are laid, that causal sets are the histories 
in a sum-over-histories quantum theory of spacetime. 

In detail, a causal set is a set C endowed with a binary 
relation -< satisfying: 

1. transitivity: ii x ^ y and y < z then x ^ z, 
Va;, y,z £ C; 

2. reflexivity: x ^ x, ^x C; 

3. acyclicity: if x ^ y and y ^ x, then x — y, Vx, y G 
C; 

4. local finiteness: Vx, z € C the set {y \ x < y ^ z} 
of elements is finite. 

Our observed continuum Lorentzian manifold, it is 
assumed, arises as an approximation to an underlying 
causal set. The partial order gives rise to the causal or- 
dering of events in the approximating continuum space- 
time, and the number of elements comprising a space- 
time region gives the volume of that region in funda- 
mental units which we take to be of order the Planck 
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volume. The above rules of correspondence give an es- 
sentially unique way to associate a class of causal sets 
to a given continuum spacetime via a process known as 
"sprinkling" defined as follows. Given a Lorcntzian man- 
ifold, {M,g), points are selected from M randomly via a 
Poisson process in which the probability measure is equal 
to the spacetime volume measure in some fundamental 
units. The selected points are the elements of a causal 
set once they have been endowed with the partial order 
induced by the spacetime causal order. The number of 
points chosen from any region of the manifold will be 
approximately equal to the volume of the region (in fun- 
damental units) up to Poisson fluctuations. For more 
details on sprinklings see the reviews mentioned above, 
for a proof of the Lorentz invariance of the process see 
We thus have a straightforward way to construct 
a causal set that could be the discrete underpinning of 
a particular continuum spacetime and, consequently, a 
starting point to develop the phenomenology of discrete 
spacetime. 

An obvious place to look for consequences of causal set 
theory is in the behaviour of particles. If the underlying 
spacetime is a discrete structure rather than a contin- 
uous manifold, free particles might no longer be able to 
follow precise timelike geodesies. Intuitively, the underly- 
ing discreteness could cause the particles to 'swerve' and 
indeed a model of particle behaviour illustrating this was 
proposed in [l^. There, a classical particle is modelled 
in the simplest possible way, as a point with no internal 
structure. The causal set, C, considered is a sprinkling 
into Minkowski spacetime and a particle trajectory con- 
sists of a chain of elements where a chain is a totally 
ordered subset of C. The trajectory is constructed itera- 
tively, where the trajectory's past determines its future, 
but only a certain proper time r/ (the "forgetting time" ) 
into the past is relevant. If the particle has reached an 
element e„ with four-momentum p„, the next element 
Cn+i is chosen such that 

• e„+i is in the causal future of e„ and within a 
proper time Tf, 

• the momentum change \pn+i — Pn\ is minimised. 

Here the momentum Pn+i is defined to be proportional 
to the vector between e„ and e„+i and on the mass shell. 
Heuristically, the trajectory tries to stay as straight as 
possible at each step. Indeed, an alternative way to define 
the process is to specify that in a frame in which the 
last two elements of the trajectory, e„_i and e„, lie on 
the vertical i-axis, the next element e„+i is chosen to be 
the one, within proper time Tf to the future of e„, such 
that the vector from e„ to e„+i is as close to vertical as 
possible. 

In this simple model the discreteness of a causal set re- 
sults in random fluctuations in the momentum of a par- 
ticle. 

One could object that this model is not intrinsic to 
the causal set as it makes use of information in the con- 



tinuum manifold to define the momentum change. How- 
ever, similar models can be defined with no reference to 
the continuum. Two such models are proposed in Sec- 
tion ini One of our main claims is that, whatever the 
microscopic model of particle motion, if it is Lorentz in- 
variant and gives rise to small random fluctuations in the 
momentum of the particle then it can be approximated 
by a continuum description as a diffusion in momentum 
space. In Section Hill we support this claim by giving the 
derivation of the diffusion equation for massive particles 
introduced in [l^. We also derive the particle diffusion 
equation in a more useful cosmic time form and without 
the original assumption of spatial homogeneity. 

In Section IIVI we explore the case of massless parti- 
cles on a causal set and obtain diffusion equations for 
the momentum of massless particles in the continuum 
approximation. Bounds are placed on the constants in 
the massless particle diffusion equation in Section |V] by 
considering the effect of momentum diffusion on the spec- 
trum of the cosmic microwave background. We will use 
units in which c = h = G = 1 — which we will refer 
to as "Planck units". Fundamental units are related to 
Planck units by a, yet to be determined factor of order 
1. Boltzmann's constant is also set to one, fc^ = 1. 



II. INTRINSIC MODELS FOR MASSIVE 
PARTICLES 

As mentioned above, the original microscopic model 
in (l^ depended on information from the continuum 
Minkowski spacetime whereas a better model ought to 
be intrinsic to the causal set itself and rely only on the 
order relation. Two slightly different intrinsic models will 
be described in this section, to give an idea of the wealth 
of possibilities available. 

We first recall some causal set definitions. Let C be a 
causal set. 

• A link is an irreducible relation, i.e. a pair of dis- 
tinct elements a, b such that a ^ b and there exists 
no distinct c such that a < c ^ b. 

• A chain is a totally ordered subset of C. An n 
chain is a chain with n elements and its length is 
n — 1, the number of links. 

• A longest chain between two elements is a chain 
whose length is maximal amongst chains between 
those endpoints. There may be more than one 
longest chain between two elements. 

• On a causal set the closest approximation we have 
to a timelike geodesic between two elements is a 
longest chain. For two causal set elements a and b 
the length of a longest chain between a and b will 
be denoted d{a,b). For sprinklings into Minkowski 
spacetime, in the asymptotic limit of large dis- 
tances, d(a, b) ~ aT where T is the proper time 
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between a and b and a is a (dimension dependent) 
constant p^ . 

• There is a hnk between elements a and b iff d(a, b) = 
1. 

• A path is a chain consisting entirely of links, i.e. a 
set of elements a-< 6-< c-< d < ... such that 
d(a, 6) = 1, d{b, c) = 1, d(c, d) = 1 . . .. 



A. Model 1 

If a dynamical rule for particle motion is to be intrinsic 
to the causal set background it can no longer refer to a 
forgetting time Tf. This instead becomes a 'forgetting 
number', an integer >> 1. In Intrinsic Model 1 a 
particle trajectory is a chain, . . . e„_2 -< e„_i -< e„ . . . 
which is determined by the following (Markov of order 2) 
process: 

Given a partial particle trajectory . . . e„_i, e„ the next 
element e„+i is chosen such that 



d(e„, e„+i) = n/, 
d(e„_i, e„+i) = 2n/, 



(see Figure 1(a) I. These requirements do not guarantee 
the existence of a unique such e„+i. However there will 
almost surely be finitely many eligible elements and we 
therefore construct the trajectory by choosing an element 
uniformly at random from these. 

The particle trajectory should swerve a little, but re- 
main approximately straight so long as nj is large, since, 
in that case, the results of Brightwell and Gregory [T3| 
show that the expected position of e„+i is close to the 
hyperboloid of points proper distance Uf /a from Cn and 
to the hyperboloid of points proper distance 2n f /a from 

In this model we can consider the trajectory as con- 
sisting of just the elements . . . e„_i, e„, e„-|-i . . . or of 
the "filled in chain" consisting of a (randomly chosen) 
longest chain (of length Uf) between e„_i and e„, an- 
other longest chain between e„ and e„+i (also length 
n/) and so on. By imposing (i(e„_i, e„+i) — 2nf we have 
forced the chain of length 2nf that we have between e„_i 
and e„+i also to be a longest chain. The trajectory is thus 
approximately geodesic over all {e„_i : e„+i} segments. 
The trajectory consisting of longest chains between e„_i 
and e„, e„ and e„+i, and e„+i and e„+2 is not, however, 
necessarily a longest chain between e„_i and e„+2- 

Possible variations on this model include choosing, at 
random, the forgetting number at each step so that the 
mean is with some fixed variance. 



B. Model 2 

The trajectory is explicitly constructed as a path in 
this model, i.e. d(e„, e„+i) = 1 for all n. Given a partial 



particle trajectory . . . e 
e„ is chosen such that 

• c^(e„-i, e„) = 1, 



n—n f 1 • 



the next element 



d(e„ 



imised, 



e«) + - ■ ■ + d{e ) is min- 



(see Figure 1(b)). Note that this minimisation docs not 



necessarily yield a unique e„, in which case we construct 
the trajectory by choosing an element uniformly at ran- 
dom from those eligible. Also, if the past trajectory has 
length less than Uf the minimisation is done over all el- 
ements available. 

Each element is linked to the previous, i.e. 
(i(e„_i,e„) = 1 so we know there exists a chain (our 
trajectory) of length nf between e„_„^. and e„. The 
maximal chain length, d{en-nf,Sn), must therefore be 
greater than or equal to n^. If we choose e„ to minimise 



d{er. 



, e„) we ask that the trajectory be as close as 



possible to geodesic between en-nj and e„ while fulfill- 
ing d(e„_i,e„) = 1. Minimising the sum of the partial 
lengths distributes the geodesic property along the path. 



III. THE CONTINUUM APPROXIMATION 
FOR MASSIVE PARTICLES 



The models described above are intrinsic to a causal set 
and as such could be used to define particle motion on any 
causal set whatsoever. For phenomenology however the 
models are only of interest when defined on a causal set 
that arises by sprinkling into four-dimensional Minkowski 
spacetime, M*, or some other physical spacetime such as 
a Friedmann-Robertson- Walker (FRW) cosmology. It is 
the central conjecture, or Hauptvermutung, of causal set 
theory that a causal set that arises by sprinkling into 
M'', is well- approximated by M'* and that the associated 
embedding of the causal set in gives us a way to derive 
spacetime physics from physics on the causal set. The 
Hauptvermutung remains to be proved - though we have 
a quite a bit of evidence for it especially in the case of flat 
spacetime (see, for example [l3lll4|) - and in the context 
of the present work it is the central phenomenological 
assumption. 

With this assumption, the models of Section |TT] defined 
on a causal set which arises from sprinkling into are 
also models of particle motion in M^. To be completely 
explicit: in both models the particle's motion on the 
causal set - a chain - defines a piecewise linear trajectory 
in M"* in the following way. By assumption, the causal 
set arises from sprinkling and so is embedded in M** (in 
such a way that the order on the elements respects the 
causal order on the embedded positions and the number 
of elements in any region is approximately the spacetime 
volume in fundamental units). Between two causal set el- 
ements in the chain of the particle's motion on the causal 
set, its spacetime trajectory through M** is a straight line 
and its four-momentum is well-defined by requiring it to 
be on the mass shell, hyperbolic space H3. 
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d{n - 1,11 + 2) > 3?!/ 
d{n + 1, ji + 2) ^ 77/ 

d(n.n + 1) — -rif 



d{n - l,n. + 1) ^ 277; 



£i(7l,77 + 2) = 277; 



ti(71 - 1,77 + 2) > 3 
£,,+2 

d(77 + 1,71 + 2) = 1 



rf(71, 71. + 1) = 1 



d(77 — 1, 77) — 77; 



d{n - 1, 77) = 1 



(b) 



FIG. 1: A trajectory constructed using (a) model 1 and (b) model 2. 



The original swerves model can be similarly converted 
into a model of piecewise linear particle worldlines in 
and we can consider each of these spacetime models as 
defining an evolution on a manifold of states. The mani- 



fold of states for the particle is . 



I3 - its position m 



space and its momentum - and the time parameter for 
the evolution is proper time. This evolution is effectively 
stochastic because although knowledge of the causal set 
makes the trajectory deterministic, we treat the causal 
set as unknown, the analogue of the state of the water 
molecules causing Brownian motion. 

If we define rimacro to be the scale of macroscopic 
physics measured in Planck units then we can demand 
a separation of scales so that 1 << n/ << nmacro- In- 
deed, the Planck scale is so small that even if we want to 
choose the scale of Large Hadron Collider physics as our 
"macro" scale (rimacro = 10^^) there's still plenty of room 
to chose Hf. The change in spacetime position which is 
bounded by nj Planck units will then be small at each 
step as will the change in momentum. 

The dynamics is therefore a stochastic evolution on 
a manifold of states and such systems are dealt with in 
general in the formalism developed by Sorkin jl^ . At the 
macroscopic scale of many (nmacro /'t-/) steps the process 
can be described approximately as a diffusion. 

Although the models described above cannot be con- 
sidered completely realistic (for example the particles are 
classical and zero-size) we claim that provided the pro- 
cess is Lorentz- and translation-invariant, it will always 
give rise to the same diffusion equation, namely the equa- 
tion written down in (l^ . As promised there, we present 
below the full derivation of this equation, supporting our 
claim that the continuum model is universal and indepen- 
dent of the discrete microscopic details. We derive the 
equation initially with the particle's proper time play- 
ing the role of independent variable; we then obtain the 
equivalent equation in terms of cosmic time, by express- 
ing both in terms of a conserved current in a certain 
space of eight dimensions. This Lorentz invariant pro- 
cess was first considered by Dudley [l^, [l3| , though one 
of his diffusion equations conflicts with ours. Without 
any imposition of Lorentz invariance, a general formal- 



ism for describing diffusion in Minkowski space was set 
up by Schay 



A. The diffusion equation for a massive particle 



We use the general formalism of [15| , which deals with 
stochastic evolution on a manifold of states. The state 
space, Ai, of the swerving particle of mass m is = 
M** X H'^, where is the mass shell. The coordinates 
on M'' are the usual Cartesians {x^}, /z = 0,1,2,3 and 
indices are raised and lowered with ry^j^, the Minkowski 
metric. The spatial coordinates on M will be written as 
{a;*}. Cartesian coordinates in momentum space are 
and whenever they are used it will be understood that 
lies on the mass shell which is the hyperboloid in mo- 



mentum space defined by Pfj,p^ + m? 
energy and p = 



0. p^ = E is the 



-Pi 



p| is the norm of the three 
^ will be written 



momentum. The three coordinates on 
abstractly as p°. We denote the coordinates on A4 col- 
lectively as — {a;^,p°} and in what follows capital 
letters A, B will be used to indicate general indices on 
M.] fi, V are indices on M^; i, j are spatial indices on M"*; 
a, h are indices on IHI'^. 

The metric on M is the product of the Minkowski met- 
ric r]^i, on and the Lobachevski metric gab on H'^. 
This is the unique Poincare invariant metric (up to an 
overall constant). The "density of states", n, plays a role 
in the formalism of [15], and by symmetry, it must be 
proportional to the volume measure on A1, so n oc ^ 
where g — det(gah). The "entropy scalar", s, is given by 
s = ln(n) (Boltzmann's constant has been set to l).[23j 

A process that undergoes stochastic evolution on a 
manifold of states, M., in time parameter T, can be de- 
scribed by a current, and a continuity equation [l^: 



dp_ 

dT 



(1) 
(2) 



Here the probability density for the system is given by 
p — p (X"^,T), a scalar density on M. The coefficients 
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K are given by 

where < • > denotes expectation value in the process in 
which the particle starts at a definite point of M. (page 
146 of [3). K^^ is a symmetric, positive semi-definite 
matrix which transforms as the components of a tensor 
on M. The coefficients are 

yA ^ / ^\ (4) 

AT-+0+\ AT / ' ^ ' 

and do not transform as a vector on Ai , but can be com- 
bined with K and the entropy scalar s to form a true 
vector u^, 

u^^v^~ dBK^"" - K^^'dBS . (5) 

The current and continuity equations can be reex- 
pressed in terms of the true vector u^: 

To find the diffusion equation for our particle process, 
therefore, we need to determine K^^ and u^. 

Requiring the equation to be Poincare invariant is a 
very stringent condition and proves to be sufficient for 
us to determine K^^ and tf^, up to the choice of one 
constant parameter. This means that the resulting equa- 
tion is very robust and independent of the details of the 
underlying particle model so long as it is Poincare invari- 
ant. 

Consider the process referred to r, proper time along 
the worldline of the particle. Then 

Ar^o+\ 2Ar / ^ ' 

Ax^ = —p^Ar at every step of the process and so 
Kf'" = iTimAr^oP''p''Ar = 0. Since K^^ is positive 
definite, this implies that K^^"^ = and the only nonzero 
components are K"'''. The only Lorentz invariant tensor 
on is proportional to the metric, g"''' and the coeffi- 
cient is independent of a;^ by translation invariance. So 
we have 

where fc > is a constant. 
Now consider 

vf^ = hm / ^\ , (9) 

Ar^O+\ At / ^ ^ 

which, by the above is w'* — /m. The components 
of the true vector are equal to because K^^^ — 0. 
There is no Lorentz invariant vector on and so = 0: 

u"^ ^ (pi' /m,Q) . (10) 



We can now write down the proper time diffusion equa- 
tion from HI) and (g]): 

|^ = *a,(.,-'^*(^))-V9,P,. (11) 

If we define a scalar 'p — pr/ ^/g we obtain the equation 
in reference ^] : 

l^^fcVl^p- V^.P, (12) 
where is the Laplacian on H'^. 



B. Diffusion in cosmic time for massive particles 



t 

i , constant r 




constant t 



FIG. 2: Particle trajectories as flowlines in X' = MT* x x R 
(where we have suppressed two spatial dimensions). 



Given an initial distribution of particles, for instance 
from an astronomical source, the above equation is not 
very useful for predicting the results of observations. 
Even if particles all leave the source at the same time 
with the same momentum, the momentum variation in- 
duced by the swerves will result in particles arriving after 
different proper times and at different observatory times. 
The proper time that elapses along the particles' world- 
lines from source to detector is not observable. To com- 
pare the swerves model with experiment and observation 
it is necessary to describe the evolution of the distribu- 
tion in time in the rest frame of our detector, which time 
we refer to as cosmic time. 

A first step in this direction was to look at the non- 
relativistic limit of the proper time diffusion equation, 
when proper time and cosmic time are comparable. The 
nonrelativistic limit in fact proves sufficient to place very 
strong bounds on the value of the diffusion constant and 
severely limit any observable effects (see [13, ) • 

In the fully relativistic case, Dowker et al. wrote down 
the diffusion equation in terms of cosmic time for the 
special case of an initially spatially homogeneous distri- 
bution [T^ l . We will now give the derivation of the cosmic 
time evolution equation for the general case of spatially 
inhomogeneous distributions. 
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The conversion between proper time and cosmic time 
is possible because both are good time parameters along 
all possible particle worldlines, which are causal. If we vi- 
sualise our diffusion process as a collection of such world- 
lines through spacctimc and momentum space, both cos- 
mic time, t — x'^ , in our chosen frame and proper time 
T increase monotonically along each trajectory. Adding 
proper time to our state space by assuming that the par- 
ticle starts at parameter r = and cosmic time t = 0, the 
process is represented by flowlines in A^' = x H"^ x K 
(see Figure [2]) and along each flowline, both r and t are 
good time parameters. The proper time diffusion equa- 
tion we have found describes the evolution of the dis- 
tribution on constant r hypersurfaces in M' . What we 
want is to obtain the diffusion equation for evolution of 
the distribution on constant t hypersurfaces integrated 
over all proper times. 

First we put t and r on an equal footing by considering 
the larger space M' and defining a new current compo- 
nent 



(13) 



If we denote coordinates on this extended space, A4' ~ 
X R, by X" = {X^, t} then the continuity equation 
^ can be written 



5a J" = 



(14) 



Using equation ([T]) (and still treating r as our time- 
parameter) we can express the t component of the current 
in terms of J"^ (a.k.a. pr). 



(15) 



where 7 = dt/dr is the usual relativistic gamma factor. 
The remaining components of the current can now be 
written in terms of J*. The spatial components are: 



m 7 



(16) 



In the case of the p components the algebra is simpler if 
we first note that we can express ([1]) in the form (cf. ([6])) 



J"" = -K 



AB 



("^-(9) 



and so 



J''{t,x\p'',T) = -kg^^ndb 



pu 



7 n 



(17) 



= -kg'^'Vad, (:^) . (18) 

The metric g"** that appears here is the Lobachcvski met- 
ric on H'^. 



Since r is unobservable we need to integrate J over r 
and we denote the integrated current by J. Integrating 
the t component of the current over proper time from 
zero to infinity, gives us the probability density on a hy- 
persurface of constant t: 



pt = J\x\p\t) 



(19) 



The components of the new current can be written: 

J\x\p'',t) EE j rdr 

J TO7 

J* ^ p^pt_ 
m 7 TO 7 ' 

J''{x\p'',t) = / J^rfr 



(20) 



-kg^'ndb ( — 
.in 



-kg^Kdh 



(21) 



\7n/ 

If we integrate the continuity equation over r we obtain 

+ dtJ' + d,r + dar ^ . (22) 

J^|r=o is zero for alH > and J'^ tends to zero as t goes 
to infinity for finite t. So for alH > we have 

dtJ' + dJ' + daJ'' = 0, (23) 

which gives the cosmic-time diffusion equation 



dpt 
dt 



-d,pt + kda[ g'^'Vadb , ^ 

m7 V \lV9 



(24) 



This is a powerful phenomenological model because it 
depends on only one parameter, the diffusion constant k. 
Data can therefore strongly constrain k. 

We note that this solves a problem posed by Dudley 
p^ . We also point out that Dudley's equation for the 
spatially homogeneous distribution on page 267 of [lB| 
is inconsistent with our equation (Equation (3.60) 

of also differs from the 1-1-1 dimensional analog of 

(El.) 



IV. MASSLESS PARTICLES 

If an underlying spacetime discreteness results in diffu- 
sion in momentum and spacetime for massive particles, 
it is interesting to consider whether a similar diffusion 
occurs for massless particles. For massive particles the 
concrete models for particle dynamics on a causal set de- 
scribed above motivated the derivation of the diffusion 
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equation (jlip . The case of massless particles on a causal 
set background is rather different. If we consider a sprin- 
kling into Minkowski spacetime, for any given element, 
p, there will almost surely be no element sprinkled on 
the future light cone of p. The analogue of the future 
light cone of p in a sprinkling of Minkowski spacetime is 
the set of all elements preceded by and linked to p. The 
elements are distributed, roughly, near the hyperboloid 
one Planck unit of proper time to the future of p. Al- 
though the whole light cone thus has a good causal set 
analogue [2J|, the easiest analogue of a null ray is only 
a single link, making it hard to see how to construct a 
discrete Markovian process that would result in a close- 
to- null trajectory. 

Modeling the propagation of massless point particles 
on a causal set as an approximately local process is there- 
fore problematic. It is hoped that in the future, the study 
of massless fields on a causal set will enable us to model 
massless particle propagation as wave packets, say. In 
the meantime, however, lack of knowledge of the exact 
nature of massless particle propagation on the discrete 
level, does not mean we cannot derive a diffusion equa- 
tion to describe the potential effect of discreteness on 
photons in the continuum approximation. We can arrive 
at a massless diffusion equation in two ways: using the 
stochastic evolution on a manifold of states procedure as 
for the massive particle case, or simply taking a to — > 
limit of the diffusion equation for massive particles. It 
turns out that the second method gives an incomplete 
result. 

The state space in the massless case differs from the 
massive case. For massive particles we had a probabil- 
ity distribution on M'* X H^. For massless particles H'' 
becomes the light cone in momentum space defined by 
PfiP^ = 0. This cone will be denoted Hg. If we assume 
that the photons under consideration are well described 
in a geometrical optics approximation so they have defi- 
nite spacetime worldlines and momenta, our state space 
will be [1^ X Hq. Since proper time vanishes along 
a lightlike worldline, it is no longer a suitable time pa- 
rameter for our diffusion process. We define, instead, an 
affine time. A, along any photon worldline by 

dxf" =p>'d\ . 

Notice that the normalization of this affine parameter is 
not arbitrary. It is fixed by its relation to the particle's 
four-momentum, or geometrically, to its de Broglie wave- 
length. Under the latter interpretation, the affine param- 
eter along a photon worldline 7 measures the area swept 
out in spacetime by a vector connecting 7 to a neighbor- 
ing null geodesic that trails it by one wavelength. 

In the massive particle case we equated the density of 
microstates, n, to the determinant of the metric on our 
state space: n cx In the massless case, the met- 

ric induced on Hg degenerates, but Hq still possesses a 
Lorentz invariant measure of volume (unique up to a con- 
stant factor). The four dimensional volume element d'^p 
of momentum space, together with the masslessness con- 



straint, p^p^ = 0, lets us construct on Hg the invariant 
volume element d'^p5{p^^p^) — d^p/2p^, i.e. n oc l/p° in 
Cartesian coordinates. It will be more useful, however, to 
work in polar coordinates on Hg: {p, 6, 0} where p is the 
magnitude of the three momentum and 6 and <j) are the 
usual polar angles in momentum space. In these coordi- 
nates, the density of states is n cx psvaO. There is also a 
(unique up to a constant factor) invariant vector field on 
Hq which is the momentum itself, i.e. the vector with 
components (p, 0,0) in polar coordinates. This is absent 
in the massive case, where the momentum vector does 
not lie in the mass shell. Finally, although there is no 
invariant metric on Hg, there is an invariant symmetric 
2-tensor, p'^p^ (unique up to a constant factor). 

We first consider the process in affine time, A. As 
with the massive case, we begin with the current and 
continuity equations, ([1]) and and determine K^^ 
and . Using the formulae ([3|) and Q with T = A we 
find 



lim p'^p'^AA 

AA^O+ 



0. 



(25) 



K is positive semidefinite so K^^'^ = 0, and finally A'"'' 
must be Lorentz invariant and translation invariant so 





/cip^p'' 



(26) 



where fci > is a constant. 

To determine we again look individually at the 
components in spacetime and momentum space. As be- 
fore, the spacetime component = hy and 
= p^ by In contrast to the massive case, there 
can be nonzero components of in the momentum 
space directions because the momentum itself is an in- 
variant vector. The momentum direction components 
are thus given by u° = ^2^", where /c2 is a constant. 
Working in polar coordinates the "position" vector on 
the cone Hg is simply (p, 0, 0) where p^ = po^. Thus 



= (p°,p^p^p^,/c2P,0,0) on 



X ml 



Substituting the forms for K^^ and into ([6]) we 
obtain the massless particle affine time equation: 



dX 



dxf^ 



k,-^ ( 9 (p. 



dE 



dE \E 



(27) 



where we have replaced p by energy E — p. 

We see that the Lorentz invariance means that any 
diffusion in photon momentum cannot change the direc- 
tion of the photon and so it always propagates on the 
light cone, at the speed of light. However the energy of 
the photon does undergo a diffusion. Notice also that 
there are two parameters, making this a less powerful 
phenomenological model than the massive particle model 
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which has a single parameter. There is not only a dif- 
fusion term but an independent drift term, arising from 
the existence of an invariant vector on Hg, and we will 
see that this leads to the existence of power law equilib- 
rium solutions. Note that taking the m ^ limit of pT|) 
would have resulted in ([27|l with fc2 = because there is 
no invariant vector in the massive case. (As is familiar 
in another context, the case of zero photon mass is thus, 
here also, a sort of singular limit of the massive case.) 



A. Cosmic time process 

Again, in order to make contact with observations, we 
need to obtain the cosmic time diffusion equation for 
massless particles, for which we use the same argument 
as in the massive case. First we assume that A = at 
t = 0. Let 



(28) 



We can then express the t component of the current J 
in terms of J^, and the remaining components of the 
current in terms of J*. 



scalar density J*, which we rename pt'. 
ot 

= -^d.,p,^{h+k,)^ + hE^, (36) 

where E ^ p is the energy. 

So, for a massless particle in a geometric optics ap- 
proximation, we expect that an underlying discreteness 
can induce fluctuations in the energy of the particle, but 
without affecting the direction of propagation. The dif- 
fusion governed by ki causes a distribution of energies 
that is initially sharply peaked to spread over time. The 
second constant k2 results in an independent drift of the 
spectrum to higher or lower energies depending on its 
sign. 

It is interesting that negative values of fc2 allow for 
power law equilibrium solutions of ([55]) . Set d^pt = 0. 
Then the equilibrium distributions satisfy 

- (fci + fc2) + k,E^ = . (37) 
This has a power law solution 



r (t,x\p'',\) = ^J* ; 
\ / p 

r{t,x\p\\) - -fcipVpsin^ab j^- 

7* 

P 



smf 



(29) 
(30) 



(31) 



Pt oc E' 



(38) 



J° is proportional to p°' and in polar coordinates the vec- 
tor = (p, 0, 0), so there is only one nonzero component 
of in the radial (energy) direction: 

Jf (i,x',p",A) = -kip^ + {2ki + k2)J* . (32) 

The afSne time of flight is unobservable so we integrate 
over it. Defining 



rit.x'^p") = / J'(^,a;^p^A)dA, (33) 
Jo 

we integrate the other current components over A to ob- 
tain 



J'(^,a;^p°) = -J* 
P 



(34) 



JP{t,x\p^) = -fcip^ + (2A:i + ^2) J*. (35) 
ap 

Imposing the continuity equation, gives us the massless 
particle cosmic time diffusion equation in terms of the 



When the parameters are such that the exponent is less 
than —2 (and so ^2 must be negative because fci is pos- 
itive) then this solution is normalisable if it is cut off at 
small energies. We conjecture that if {2ki -\- fc2)/fci < —2 
any normalised distribution will tend at late times to this 
power law equilibrium solution at large energies. This 
is interesting because physical processes that result in 
power law distributions across a wide energy range are 
few and far between — the Fermi mechanism of statisti- 
cal acceleration of charged particles by random magnetic 
fields, proposed as the source of high energy cosmic rays, 
is the only well-known mechanism. 

Placing bounds on the parameters fci and /c2 is the next 
step. 



V. BOUNDING THE CONSTANTS fei AND k2 

In developing a phenomenological model, one aims to 
provide a model for currently unexplained observations 
or suggest new observations that might be made to test 
a theory. But before proposing new observations, one 
should of course constrain one's model as tightly as pos- 
sible, based on what is already known. Our model has 
two parameters: a positive diffusion constant ki and a 
"drift" constant k2 , which may be either positive or neg- 
ative. To place the strongest bounds on the values of 
these parameters, it seems sensible to look at photons 
that have been travelling for a very long time and thus 
have had the maximum possibility to experience any un- 
derlying discreteness. The cosmic microwave background 
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(CMB) seems an ideal testing ground in this sense. Not 
only are its photons the "oldest" we can observe, but 
its spectrum has been determined with great precision. 
Most of the photons in the CMB have been "free stream- 
ing" for approximately 13.7 billion years, or on the order 
of 10^" Planck times. When the universe became trans- 
parent at recombination, they would have had a black- 
body spectrum with a temperature 3000 K (see for ex- 
ample [20|]). Current observations yield a temperature of 
2.728 ± 0.004 iiT and measure the spectrum to be Planck- 
ian (blackbody) over the 2 — 21cto^^ frequency range 
to within a weighted rms deviation of only 50 parts per 
million (ppm) of the peak brightness [21|. Since our dif- 
fusion would have distorted the energy distribution, the 
fact that the CMB photons have travelled so far but re- 
mained so perfectly thermal will allow us to constrain our 
parameters very tightly. 



A. Simulations 

Our derivation of the massless cosmic time diffusion 
equation assumed spacetime to be Minkowskian. There- 
fore we will first consider a simplified model that ignores 
the expansion of the universe, and consequently assumes 
that, in the absence of diffusion, the temperature of the 
CMB would remain constant from the surface of last scat- 
tering to today. This will give us an order of magnitude 
bound on the parameters. In Section I VII the cosmic ex- 
pansion will be incorporated. 

The initial Planckian spectrum, expressed as a num- 
ber density of photons per unit spatial volume per unit 
energy, is 



p{E,t = 0) = Stt 



exp (I) - 1 



(39) 



with a temperature T — 2.728/^. According to our 
model, this distribution evolves via the homogeneous 
massless cosmic time diffusion equation 



^ dpt 



kiE 



dE^ 



(40) 



Using the MATLAB numerical pde solver pdepe, this equa- 
tion was integrated over a time interval equal to that 
since the surface of last scattering. 

Although only the 2— 21cto~^ region of the spectrum is 
needed to compare with the reported rms deviation, these 
evolutions were run over a larger range of frequencies to 
capture more of the spectrum and allow the implementa- 
tion of a boundary condition at E' = 0. What boundary 
condition is appropriate? What happens to a photon as 
its momentum approaches zero? Do photons leak away 
through the tip of the null cone in momentum space? 
Physically, the photon concept employed by our model 
breaks down as the wavelength tends to infinity, because 
the geometrical optics approximation fails. (Moreover, 
our affine parameter A fails to be well defined physically. 



since it reaches infinity in a finite time \fE ^ 0, and thus 
cannot remain approximately constant over the photon 
wave-packet.) This suggests that the so called "absorbing 
boundary condition", p{E) — 0, is appropriate at -E = 0, 
and this is what was used in all our simulations. In fact, 
the current is 



J = (2fci + k2)pt " kiEdpt/dE , 



(41) 



so (as long as dpt/dE remains finite) any linear combina- 
tion of pt = with the "reflecting boundary condition" , 
J = 0, is equivalent at E = 0. 

The evolved spectrum was converted from a number 
density per unit volume per unit frequency to a spectral 
radiance — energy per unit area per unit time per unit 
frequency per steradian — as used in the analysis of the 
COBE FIRAS data. This allows us to compare the devi- 
ation from Planckian with the quoted 50ppm of the peak 
brightness. 

A Planck spectrum was fit to the evolved spectral radi- 
ance using the least squares method. By looking for the 
best fit Planck spectrum rather than comparing with the 
initial 2.728i4r spectrum, we allowed for the possibility 
that the diffusion changes the temperature of the CMB 
in a way that may be reconciled with observation. As 
it happens, we found that the temperature of the best 
fit Planck spectrum was very close to the initial tem- 
perature in cases where the deviation is within the al- 
lowed tolerance. For example the choice of parameters 
ki = 5 X 10^^^ and k2 = 1 x 10^^^ gives a best fit tem- 
perature of 2.7281i4r, indistinguishable from the current 
observed temperature of 2.728 ±0.004^. Finally the rms 
deviation between the fitted Planckian spectrum and the 
evolved spectrum in the 2— 21cto~^ frequency range (en- 
ergy range 4 x 10^^'^ — 4 x 10^^^ J) was calculated with 
all points weighted equally. This result was compared 
to the allowed tolerance of 50 parts per million of the 
peak brightness. This process was repeated for a range 
of values of the parameters fci and k2. 



B. Results 

We first place bounds on the diffusion and drift con- 
stants separately, varying fci with k2 — and varying k2 
with fci = 0. 

When fci = we can solve the equation exactly: 



PtiE,t)^ po{E-k2t). 



(42) 



so the spectrum just translates at a constant speed. For 
k2 negative, this is inconsistent with the boundary condi- 
tion p — ai E — 0. However, in this case one can imple- 
ment an absorbing boundary condition trivially: simply 
cut off the translated distribution at E = 0. This is what 



we did to generate the solution plotted in Figure 3(b) 



One might be concerned that deviations within the al- 
lowed tolerance would be so small as to approach the 
level of the numerical errors in the simulations. The 
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FIG. 3: The rms deviation between the evolved spectrum and a bestfit Planck spectrum (as a proportion of the spectrum 
peak): (a) varying fei with fc2 = 0, and (b) varying k2 with fci = 0. 




FIG. 4: The rms deviation of the simulated spectrum from 
Planckian as a proportion of the spectrum peak, varying both 
ki and k2. Values of ki and k2 within the 5e— 05 contour give 
a spectrum that is Planckian to within 50ppm of the peak. 



exact solution for ki — provides us with a means of 
demonstrating that this is not the case [l^. When wc 
compare the exact solution with the numerical solution 
for fci = the errors introduced by the numerical inte- 
gration can be seen to be several orders of magnitude 
smaller than the deviation from Planckian. For exam- 
ple if fc2 = 4 X 10~^^ the rms deviation from the best 
fit Planck spectrum is 5 x 10"^"^ (5 x 10~^ peak bright- 
ness) for both the exact and the numerical solution. The 
rms deviation between the exact and numerical solution 
is 4 X 10~^°^. If fc2 = — 4 X 10^^^ the rms deviation from 
the best fit Planck spectrum is also 5 x 10^^"^ (5 x 10"'^ 
peak brightness) while the deviation between the exact 
and numerical solutions is again 4 x 10"^°''. This also 
demonstrates that the p = boundary condition we im- 
posed on the numerical solution, although inconsistent 



with the exact solution when k2 < 0, does not intro- 
duce noticeable errors for the values of ki that we are 
concerned with. 

When k2 — with /ci > we can only solve the equa- 
tion numerically. The results for both cases are displayed 
in Figures [3(a)] and [3(b)| We see that the deviation from 
Planckian increases approximately linearly with increas- 
ing magnitude of the parameters. (Notice that figure [3(b)] 
was drawn from the exact solution, the graph taken from 
the numerical solution is indistinguishable.) The simu- 
lations suggest that for the deviation from Planckian of 
the CMB to be within the allowed 5 x 10~^ of the peak 
brightness the diffusion constant ki must be less than ap- 
proximately 7 X 10"^"^ if k2 = 0, and the drift parameter 
/c2 must fall within the range — 4x 10~^^ < k2 < 4x 10 
if ki ~ 0. Converting to SI units we have the bounds 
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ki <3x IQ-*"^ kgrn'^s^^ , 



(43) 

-1 X 10"^^ < A:2 < 1 X lO^^^/cgm^s^^ (44) 

Similar bounds apply when we let both ki and k2 be 
nonzero. The general situation is displayed in Figure IH 
from which one can read off the values of fci and ^2 for 
which the deviation from blackbody is less than 5 x 10~^ 
of the peak brightness when we allow both constants to 
vary. 

In the units used here, the bounds on the parameters 
are very small. However we can get a handle on where 
these numbers come from by rescaling the energy, setting 
E' = sE with s chosen so that sT = 1 when T is the CMB 
temperature. This means that s ~ 10"^^ in Planck units. 
We rescale p' = p/s so the initial spectrum is: 



1 



E' 



,1 



(45) 



where T' = sT = 1. If fc2 



0, then we can also rescale 
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the time, setting t' = skit to obtain the diffusion equa- 
tion 



dp' 



E' 



(46) 



If we now evolve p' until it differs from p'q by 50 ppm 
and take the value, t'j: of t' when this happens, t'j: must, 
for consistency with the data, be greater than or equal to 
s kit where t is the age of the universe, and so in Planck 
units fci < 10~^'^10~^^t'j. We see that the order of magni- 
tude bound found above will result if t'j: ^ lO"**, which is 
indeed about the (rescaled) time at which one would have 
expected the deviation to reach 50 ppm. A similar order 
of magnitude estimate follows from the geometric inter- 
pretation of our afhne parameter A as an area, if one notes 
that the product of the photon wavelength 1cm) with 
the Hubble radius is around lO^^lO^o - 10^^ pianck 
units. 



VI. EXPANDING UNIVERSE 

In Section |V] we ignored the effect of the expansion of 
the cosmos on the CMB and assumed that it remained 
at a temperature of ^ 2.7 K from the surface of last scat- 
tering to today. This is of course not the case. At the 
surface of last scattering the CMB had a temperature of 
about 3000K. As the universe expanded the individual 
photons were stretched along with the space, and cor- 
respondingly diluted, leaving us with the 2.7K spectrum 
observed today. We will now show that the expansion has 
essentially no effect on our model in the sense that the 
distribution in the expanding universe can be deduced 
easily from the nonexpanding one and that the bounds 
derived from the nonexpanding simulation change only 
slightly. 

The redshifting effect of the expansion (but not the 
dilution) can be added to the model by adding to u a 
vector which has a single component in the E direction: 



. p dE „d 
Av^ = — = -E- 
dt a 



(47) 



where a{t) is the cosmic scale factor. This changes the 
continuity equation ^ to 



at 



(48) 



-^d,pt -{ki+k2)-^ 



A solution of this equation, for fci = A:2 = is 

a3 ^2 



po {E, t) = 8^ 



3 -E a 
ag e To ao _ 1 



(50) 



where ap is the scale factor at time t^. If we multiply this 

distribution by which dilutes the photons according 
to the expansion, it becomes exactly the Planck distribu- 
tion for temperature T = Tq ^ . 

If we define a new variable E = -^E and a new density 

function p{E) = ^p{E) (this being just the transforma- 
tion of a scalar density under a rescaling of coordinates: 
pdE = pdE) the distribution po{E,t) ((50)) becomes 



= 87r- 



E^ 



exp 



(*) 



1 



(51) 



which is constant in time. 

We now transform our diffusion equation to the 
rescaled quantities p and E. 

Starting with (|49|) . we have: 



T rrr, I ^ ClE 8 \ , ^ 

LHS = Ut + ap) 

\dt a dE 



(52) 



= ap + ap + dEp' (53) 
RHS = -{ki + k2)a {ap}' + kio^Ep" + d (pE^ (54) 
= -(fci + k2)a^p' + kio^Ep" + dp + dEp' , (55) 

where dot denotes time derivative and prime denotes 
derivative with respect to E. This gives 



fi, , 1 \ ^P , I, ^^^P 
— = -(fci + k2)a—^ + kiaE^^. 

dt ' ' dE dE^ 



Choosing t' such that 



dt' 
dt 



a, we obtain 



dp ,, , , 9 _ , ~ 9^ _ 

- = -iki+k2)^p+kiE^^p, 



(56) 



(57) 



which is the same as (j40p . the nonexpanding diffusion 
equation. 

That we can find expanding solutions from static ones 
is due to the scale-invariance of the null cone Hq : its 
geometrical structures are invariant under E ^ E ^ 
const X E. 

For a matter dominated FRW universe a ~ t'^^^ i.e. 
a{t) — t^/'^/tg^"^, where to is the current value of t (and 
the current value of a is 1). We have ^ = a which 
integrates to 



t ~ -— + const 



(58) 



If the range for t is 10™ then the range for is 3/5 of this. 
So the simulations we would need to do for the expanding 
case are the same as for the nonexpanding case but for 
only 3/5 of the time. This doesn't affect the order of 
magnitude of the bounds. 
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VII. DISCUSSION 

The work presented here illustrates the familiar fact 
that considerations of symmetry can bring forth a fairly 
unique phenomenological model, even when relatively lit- 
tle is known about the deeper reality the model is meant 
to represent. Starting from the assumption of an un- 
derlying spatiotemporal discreteness that nevertheless re- 
spects Lorentz invariance in the continuum approxima- 
tion, we argued that particle momenta would be subject 
to stochastic variations, and that if these variations were 
small, their effects would be describable on large scales 
as a diffusion in momentum space. The assumption of 
Lorentz symmetry lends the resulting models their power 
(by limiting the number of parameters), and it sets them 
apart from the majority of quantum gravity phenomeno- 
logical models, which break Lorentz invariance. 

For particles without internal degrees of freedom, we 
have seen that even in the absence of a definite micro- 
scopic theory, an effective diffusion model can be derived 
based on the assumed invariance alone. One can also 
imagine applying this idea more generally, including for 
example the polarisation of photons or neutrinos. 

In the case of massive particles, if one of the ex- 
plicit microscopic models is fixed upon, then the diffu- 
sion strength, k, will be a function of the forgetting time 
(number). This forgetting time sets the scale shorter 
than which the dynamics is nonlocal: at much larger 
scales the model is effectively local. In more realistic, 
more quantal models, the diffusion scale might also de- 
pend on such dimensionless numbers as the ratio of the 
mass of the particle to the Planck mass and properties of 
the particle's wave packet. The same possibilities exist 
for the massless case. Thus we would expect the diffusion 
and drift parameters, fci and k2, to depend on some non- 
locality scale in the underlying physics, and they could 
also depend on features of the wave packet associated 
with the photon, for example the ratio of the (peak) 
wavelength to the length or the packet. In seeking an 
underlying model of photons, the Lorentz invariant, non- 



local D'Alembertian that has recently been discovered for 
scalar field propagation on causal set backgrounds [1, 
could be valuable. Using it to evolve a wave packet of a 
massless scalar field, one could ask whether the resulting 
propagation exhibited any momentum diffusion or drift, 
and if so, what sets the scale of these phenomena. 

The parameters of our model are constrained by the 
blackbody character of the CMB radiation. Since most 
observational astrophysics and cosmology relics on elec- 
tromagnetic radiation, there are a host of other observa- 
tions that could also be brought to bear, given that our 
model entails a broadening of spectral lines as well as a 
distance-dependent shift in energy. For example, if the 
diffusion constant were set to zero, it would be easy to 
work out how the drift would affect absorption spectra 
from distant objects. It seems likely however that the 
bounds set here will be among the most stringent. 

The models discussed here describe free point particles. 
Although we expect that composite objects would be less 
affected by the underlying discreteness (for example, a 
helium ion would swerve less than a proton) we can not 
make conclusive statements without a causal set model 
for interacting particles. 
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